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Motivated by recent experiments on BaCuSi206, we investigate magnetic excitations and quan- 
tum phase transitions, driven either by pressure or magnetic field, of layered dimer magnets with 
inter-layer frustration. We consider two scenarios, (A) a lattice with one dimer per unit cell and 
perfect inter-layer frustration, and (B) an enlarged unit cell with inequivalent layers, with and 
without perfect frustration. In all situations, the critical behavior at asymptotically low tempera- 
tures is three-dimensional, but the corresponding crossover scale may be tiny. Magnetic ordering 
in case (B) can be discussed in terms of two condensates; remarkably, perfect frustration renders 
the proximity effect ineffective. Then, the ordering transition will be generically split, with clear 
signatures in measurable properties. Using a generalized bond-operator method, we calculate the 
low-temperature magnetic properties in the paramagnetic and antiferromagnetic phases. Based on 
the available experimental data on BaCuSi2 06, we propose that scenario (B) with inequivalent lay- 
ers and imperfect frustration is realized in this material, likely with an additional modulation of the 
inter-layer coupling along the c axis. 



I. INTRODUCTION 

Magnetic quantum phase transitions (QPT) constitute 
an intense area of current condensed matter research. 
Similar to classical phase transitions, the universal criti- 
cal properties of continuous QPT are determined by the 
symmetry of the order parameter and the number of 
space dimensions, d, but in addition the order-parameter 
dynamics plays a distinctive role. The criterion for mean- 
field behavior, being d > A for classical magnetic transi- 
tions, changes to d-\-z > 4 where z is the dynamical expo- 
nent. Thus, realizing quantum critical behavior beyond 
mean-field often requires low-dimensional systems, with 
effectively one-dimensional (Id) or two-dimensional (2d) 
behavior. Such reduced dimensionality may be achieved 
in structured materials, consisting of chains or planes 
with a weak three-dimensional (3d) coupling: then, there 
exists an energy scale Ez where the behavior crosses over 
from 3d at low energies to Id or 2d at higher energies. 

A number of recent experiments suggested a differ- 
ent route to reduced dimensionality near magnetic QPT, 
namely through geometric frustration. One group of 
experiments are those on heavy-fermion metals, like 
CeCug-icAua; and YbRh2Si2, undergoing a transition to- 
wards an antiferromagnetic metallic state, with proper- 
ties inconsistent with the standard theory of 3d mag- 
netic QPT in metals. ^"'^ Some of the proposed theoretical 
explanations®'^" are based on the assumption of the spin 
fluctuations being two-dimensional. Remarkably, indica- 
tions for 2d critical fluctuations have indeed been found 
in neutron scattering® on CeCug-ajAu:! at the critical 
concentration Xc — 0.1. 

More recently, experiments on Mott-insulating quan- 
tum paramagnets consisting of coupled dimers of spins 
1/2 have been performed. Magnetic QPT can be 
driven either by application of a magnetic field or of 
pressure. ^^"^® The field-driven zero-temperature transi- 
tion, from a paramagnet to an antiferromagnet with 



XY order at a field H = Hd, belongs to the uni- 
versality class of the dilute Bose gas, and the corre- 
sponding finite-temperature transition has been termed 
as "Bose-Einstein condensation of magnons" . In the ma- 
terial BaCuSi206,^^'^^ the finite-temperature transition 
line has been found "'^^ to follow Tc oc {H — Hd)^ , with a 
shift exponent ^ = 1 characteristic of a two-dimensional 
quantum critical point (QCP), in a temperature range of 
30mK < T < IK. As BaCuSi206 consists of layers of 
Cu dimers with a fully frustrated inter-layer coupling, i.e. 
a body-centered tetragonal (bet) structure of dimers, the 
results have been interpreted as dimensional reduction at 
a QCP arising from geometric frustration. ^®'^^ Wc have 
recently shown^*^ that frustration, while in general not 
leading to 2d behavior at asymptotically low energies, 
can strongly suppress the 3d crossover scale Ez- More 
precisely, Ez is reduced to Ez cx Jz/J^, compared to an 
unfrustrated situation with Ez (x Jz] here J and Jz are 
characteristic in-plane and inter-plane coupling scales. 
(We note that, at very low energies, additional ingre- 
dients like magnetic anisotropies and coupling to nuclear 
spins will modify the critical behavior of a real material.) 

Returning to the compound BaCuSi206, detailed neu- 
tron scattering experiments^^ have revealed the existence 
of at least two magnetic "triplon" modes in zero field, in- 
dicating the presence of an enlarged unit cell with in- 
equivalent dimers. In fact, the material undergoes a 
structural phase transition at a temperature of 100 K, be- 
low which the magnetic layers are no longer equivalent. 
In such a situation, the application of a field will close the 
gap of the lowest triplon mode at Hd , but the conden- 
sate will be strongly inhomogeneous along the c axis, and 
hence effectively two-dimensional. An inhomogeneous 
condensate for fields above Hd has indeed been found in 
a recent NMR experiment^'^ on BaCuSi206 which shows 
the existence of (at least) two inequivalent Cu sites with 
distinct magnetizations. 

This discussion suggests two distinct routes towards 
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FIG. 1: a) bet lattice structure of dimers, with couplings J 
(intra-dimer) , J' (intra-layer) , and Jz (inter- layer). b,c) Illus- 
trations of possible symmetry breaking. In b) the inter-layer 
couplings are modulated within a unit cell such that the 
two diagonals are inequivalent. In c) the inter-layer couplings 
are further modulated along the c axis, enlarging the unit cell. 

critical points with reduced dimensionality in layered 
quantum magnets: (A) equivalent layers with perfect 
inter-layer frustration - this is the scenario proposed in 
Ref. 18 and investigated theoretically in Refs. 19,20,24; 
and (B) inequivalent layers^^^^""^ - in this situation perfect 
frustration may still be present, but is not a required in- 
gredient. 

The purpose of this paper is twofold: First, we give a 
detailed discussion of situation (A), augmenting the re- 
port given in Ref. 20. Second, we study a model with two 
sets of inequivalent layers, i.e., two dimers per unit cell, 
aiming at a semi-quantitative understanding of situation 
(B). We shall show that qualitative differences between 
the cases with and without perfect frustration arise, al- 
lowing for a clear-cut experimental distinction. As in 
Ref. 20, our general arguments are based on an anal- 
ysis of order-parameter field theories, and quantitative 
calculations are performed using variants of the bond- 
operator approach, in both the paramagnetic and anti- 
ferromagnetic phases. We shall be primarily interested in 
the mechanisms which suppress the 3d crossover scale to 
small values. The physics below this scale shall not be ad- 
dressed in detail, as e.g. material-dependent anisotropics 
(not included in our analysis) will become important. 



A. Main results 

In the following, we summarize our main results which 
strictly apply to a bet system of dimers (Fig. 1), but 
apply with minor modifications in general to magnetic 
systems with frustrated inter-layer interaction: 

(i) At the magnetic quantum critical points, the 
asymptotic low-energy physics below a scale is in gen- 
eral three-dimensional. However, the crossover scale Ez 
is strongly reduced compared to unfrustrated layered sys- 
tems, due to either (A) frustration, or (B) a layer "mis- 
match" arising from inequivalent layers. (An exception is 
the high-field quantum phase transition between a canted 
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FIG. 2: Zero-temperature magnetization vs. applied field in 
the coupled-dimer model near Hci at T = 0, calculated using 
bond operators (Sec. IV), for three scenarios: (A) equivalent 
layers and perfect frustration (dashed), (Bl) inequivalent lay- 
ers and perfect frustration (dash-dot), (B2) inequivalent lay- 
ers and imperfect frustration (solid). Case (Bl) clearly shows 
the secondary transition as a kink in the magnetization curve, 
whereas this feature is smeared out due to the proximity ef- 
fect in case (B2). Details will be discussed in Sec. VII, the 
parameter values are those in Figs. 9,16,17 below. We caution 
the reader that logarithmic corrections, not captured by our 
approximation, may somewhat modify the shape of the mag- 
netization curve close to the critical field, but the existence of 
a sharp kink in (Bl) will be unaffected. 

and a fully polarized state at field i?c2, where perfect frus- 
tration renders the low-energy physics two-dimensional, 
due to the lack of scattering processes at T = 0.^^) 

(ii) Perfect frustration implies the absence of a linear 
coupling between two distinct magnetic condensates (on 
the even and odd layers), leading to an additional Z2 
symmetry which is spontaneoulsy broken in the ordered 
phases. For inequivalent layers, the two magnetic con- 
densates will be established at different critical points, 
i.e., the magnetic transition is generically split due to 
the absence of a proximity effect. This feature persists 
to finite temperatures and allows a distinction between 
perfect and imperfect frustration - in the latter case, the 
secondary transition is smeared out. 

(iii) Applied to the field-driven QPT in BaCuSi206, 
our results suggest that scenario (B), with inequivalent 
layers and incomplete frustration, is closest to the exper- 
imental situation. To demonstrate the qualitative differ- 
ences, we plot in Fig. 2 the magnetization curves close to 
Hci for three prominent situations (equivalent layers with 
perfect frustration, inequivalent layers with and without 
perfect frustration). Finally, to reconcile the rather large 
magnetic proximity effect as evidenced by NMR^'^ with 
the tiny 3d crossover scale^^ we propose a vertically mod- 
ulated inter-layer coupling which leads to a further sup- 
pression of 3d behavior. 

We note that a recent paper by Batista et al}^ stud- 
ied in detail a hard-core boson model on a bet lattice - 
this is the appropriate model at the Hc2 transition, and 
approximately applies also in the vicinity of Hd, how- 
ever, there it does not capture the interaction-generated 
3d dispersion. 
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B. Outline 

The bulk of the paper is organized as fohows: In 
Sec. II we introduce tlic coupled-dimer model to be 
studied in this paper, and discuss ways of breaking 
the bet lattice symmetries. In Sec. Ill we develop an 
ordcr-paranicter description of the magnetically ordered 
phases, and present a comprehensive symmetry analysis. 
For quantitative calculations at zero temperature, we em- 
ploy the bond-operator method, summarized in Sec. IV. 
Sees. V and VI contain the central results of our paper, 
both for scenario (A) with equivalent layers and full frus- 
tration (Sec. V) and for scenario (B) with inequivalent 
layers (Sec. VI). We discuss the structure of the phase 
diagram, the nature of the phase transitions, various en- 
ergy scales relevant for thermodynamics as well as for the 
magnetic excitations, and zero-temperature observables 
like e.g. magnetizations. Finally, Sec. VII applies the re- 
sults to BaCuSi206, by deriving theoretical constraints 
from various experimental results: This leads us to pro- 
pose a scenario of inequivalent layers with partially frus- 
trated and vertically modulated inter-layer interaction. 
A brief outlook concludes the paper. 

Readers primarily interested in our conclusions regard- 
ing BaCuSi206 may directly jump ahead to Sec. VII (af- 
ter having glanced at the notations in Sec. II). 



II. MODEL 

In this paper, we will concentrate on bet lattices. 
Fig. 1, consisting of two interpenetrating tetragonal sub- 
systems. Assuming magnetic moments to be located on 
the sites of this lattice, with nearest-neighbor antifer- 
romagnetic couplings, the in-plane order will be Neel- 
like, but the coupling between adjacent layers is fully 
frustrated. The geometric frustration can be seen at 
the single-particle level: For a tight-binding model with 
nearest-neighbor couplings t (intra-layer) and tz (inter- 
layer), the single-particle dispersion is given by e,- = 
2t{cosqx + cosfjj,) + cos((j',,;/2) cos((7j,/2) cosfj^. For 
positive t and small tz , the minimum of the dispersion is 
at in- plane wavevector ^| = (tt, tt), where it is indepen- 
dent of cjz - this is the result of inter-layer frustration. 

With the application to coupled-dimer materials like 
BaCuSi2 06 in mind, we will consider a Heisenberg 
Hamiltonian for dimers of spins 1/2. We decompose the 
Hamiltonian into in-plane and inter-plane parts, 

n = Hi\+nz. (1) 

The in-plane part is unfrustrated and reads 

in {ij)nm 

— Sjnm (2) 

inm 



where m = 1, 2 labels the spins of each dimer, i, j are the 
dimer site indices in each layer, and n the layer index. 
J„ and are the antiferromagnetic intra-dimer and in- 
plane inter-dimer couplings, and H is an external uniform 
field. 

A rather general form of the inter-layer coupling reads 

nj \ ^ jrimm Ac 3 

f^z — / ^ 'J z ^inm * '^i-t-A,n-i-l,m' 

(3) 

inmra' 

Jz is the frustrated coupling between adjacent layers, 
whereas Jzz represents an unfrustrated coupling between 
second-neighbor layers. The runs over four spatial 
diagonals (A = 1,...,4) such that the sites [in) and 
(z + A,n + 1) are nearest neighbors in z direction, i.e., 
with distance (±a/2, ±a/2, c) where a and c are the lat- 
tice constants (which are set to unity in the following). 
Note that J^""™ ^ has some specific inter-dimer struc- 
ture given by the mm' dependence, and frustration can 
be broken due to the A dependence, see below. The 
mm' dependence requires discussion, as different phys- 
ical processes are determined by different combinations 
of J^"* : The "bare dispersion" of the triplet excitations 
is determined by J2z = + Jz^ — Jl"^ — Jz^, whereas 
the combination J4Z = Jl^ + ^ -I- J^^ -I- J^^ enters in 
the interaction vertex. If we label the lower (upper) spin 
in each dimer with 1 (2), Fig. la, the geometry suggests 
that Jz^ is dominant and antiferromagnetic. Hence, in 
what follows we shall assume — J2z = Jiz = Jz > unless 
otherwise noted. 

Occasionally we will refer to a possible biquadratic 
inter-layer exchange term: 

f- (4) 

As we have shown in Ref. 20, an effective second- 
neighbor inter-layer coupling is always generated (cx J^) 
through interaction processes, even if a "bare" Jzz is ab- 
sent in an idealized model Hamiltonian. (The same ap- 
plies to the biquadratic term JcoU, with Jcou oc — .) 
In the microscopic calculations to be presented in this 
paper, we set Jzz = JmW = unless otherwise noted. 

A. Equivalent layers 

For an ideal bet lattice structure, all layers are equiv- 
alent, Jn = J, J'n = J' , and the inter-layer coupling is 
fully frustrated, i.e., the four diagonal bonds are equal, 
jnmm A _ jmm _ Importantly, the second-neighbor cou- 
pling in z direction is allowed by symmetry,^" Jzz in 
general. 

Distortions of the ideal bet lattice may or may not 
enlarge the unit cell. In the latter case, particularly in- 
teresting are distortions which break the perfect inter- 
layer frustration. Those correspond to a A dependence of 
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jnmm A rcndcrs the four links between one site and 

its neighbors in z direction inequivalent. The simplest 
symmetry breaking leads to different couplings along the 
two diagonals, see Fig. lb. 



B. Inequivalent layers 

Various lattice distortions can occur which enlarge the 
unit cell. This is likely the case in the low-temperature 
phase of BaCuSi2 06, which, however, to our knowledge, 
has not been fully characterized to date. For simplic- 
ity and motivated by the BaCuSi206 neutron scattering 
results, we will assume that the distortions preserve 
the tendency towards commensurate in-plane ordering 
at wavevcctor (tt, tt), but we allow for layca-depcndent 
in-plane couplings J, J' and possibly broken inter-layer 
frustration. 

Importantly, within the two sets of "even" and "odd" 
layers (i.e. layers with even and odd n) a unfrustrated 
(albeit small) coupling will be present through J^z , irre- 
spective of the presence or absence of inter-layer frustra- 
tion. Thus the simplest and physically most interesting 
scenario is one where all even planes are identical, as are 
the odd ones, i.e., we have a lattice still consisting of 
two tetragonal subsystems, labeled A (n even) and B (n 
odd). Then the couplings take the form: 

•hn = Ja, J2n — J'a^ -hn+l = Jb, J2n+1 = '^B- i^) 

This also implies two different second-neighbor vertical 

COUphngS, J2« = J/^, J2n+1 = jB, 

The coupling between adjacent layers may still have 
the full symmetry, J"""" = J™™ , implying perfect 
frustration. Alternatively, the distortion can break the 
frustration. The simplest situation with inequivalent di- 
agonals. Fig. lb, can be described by couplings ac- 
cording to 

j2nmm'A _ I JTai ^ 0*^^ 

I JTJS' ^ even ' 

T2ra+l,mm'A _ \ J^Bl ^ '•"^^ (f:\ 

\ J^g A even " ^ ^ 

Here we have also allowed for a modulation of the Jz cou- 
pling along the c axis, as illustrated in Fig. Ic - this will 
turn out to be relevant for our analysis of BaCuSi206. 
(Note that Ja 7^ Jb, J'a J'b the one hand and 
JzA 7^ JzB on the other hand break different mirror sym- 
metries of the bet lattice.) 

All scenarios included above lead to a unit cell contain- 
ing no more than two dimcrs (one from the even and one 
from the odd layers). A further enlarged unit cell, i.e., 
multiple dimers within each family of layers, only lead 
to quantitative modifications of the overall physics, and 
will be discussed towards the end of the paper. 



C. Phases 

To set the stage, we sketch the phase diagram of the 
coupled-dimer model Ti, (ignoring the complications aris- 
ing from structural distortions of the ideal bet lattice). 
At zero temperature the phase diagram is similar to that 
of the well-studied bilayer Heisenberg model, see e.g. 
Fig. 1 of Ref. 27. 

For J J',\Jz\, the zero-field ground state of Ti is 
a paramagnetic singlet, with elementary gapped triplet 
excitations. The number of excitation branches is equal 
to the number of dimers per unit cell. 

If J' dominates, an antiferromagnetic phase with bro- 
ken SU(2) symmetry and in-plane ordering wavevcc- 
tor (tTjTt) is established. In the case of perfect inter- 
layer frustration, layers with classical moments would 
be uncoupled, leaving an infinite ground-state manifold. 
This degeneracy is lifted for quantum spins by a zero- 
temperature order-from-disorder mechanism, which es- 
tablishes true 3d order (Sec. VB). In contrast, for large 
\Jz\ the in-plane order is ferromagnetic, with the inter- 
plane order depending on the sign of Jz- 

Applying a field to the large-,/ quantum paramagnet 
leads to a Zeeman splitting of the triplet excitations. At 
a critical field Hd, the gap of the lowest mode closes, 
and a quantum phase transition to a gapless canted phase 
occurs. (For modifications due to inequivalent dimcrs see 
Sec. VI.) The canted phase has a broken U(l) symmetry, 
corresponding to XY order perpendicular to the field. 
Upon further increasing the field, the system is driven 
into a fully polarized state at Hc2- 



III. SYMMETRIES AND ORDER-PARAMETER 
FIELD THEORY 

In this section, we present the ingredients for a Landau- 
Ginzburg- Wilson description of the degrees of freedom of 
the coupled-dimer model H. We start with antiferromag- 
netic fluctuations on a paramagnetic background. 

A remark of caution is in order: The formal deriva- 
tion of an effective theory for the staggered magnetiza- 
tion involves integrating out the degrees of freedom as- 
sociated with the uniform magnetization, see e.g. Ref. 4. 
(The "local" staggered and uniform magnetizations cor- 
respond to Si — S2 and Si + S2, respectively, in terms 
of the two spins Si, §2 of a dimer in H (1).) However, 
there exist processes which are strictly forbidden within 
the effective theory for the staggered magnetization, but 
exist in the full theory. One example is the two-particle 
(as opposed to three-particle) decay of triplons discussed 
recently in Refs. 30,31; another example is the effective 
vertical second-neighbor hopping of order in the model 
Ti (1), to be described below in Sec. VA. 



5 



A. Antiferromagnetic order parameter 

A continuum description of antiferromagnetic fluctu- 
ations on bet lattices requires eare due to the geomet- 
ric frustration. In-plane magnetic fluctuations are strong 
and centered around wavevector Q = (7r,7r), but the ef- 
fective inter-layer coupling is frustrated and weak. Hence 
there will be a large regime of energies or temperatures 
with no well-defined order along the c axis. Therefore we 
only take the continuum limit w.r.t. the in-plane coordi- 
nates, but keep the discrete layer index n. We thus define 
an order parameter (j!)„(f|[), where r|| is the in-plane co- 
ordinate, with the local magnetization operator given by 
m„(r||) = exjpiiQ ■ r||)0„(r||). As outlined above, we will 
allow for unequal layers, but assume a two-dimer unit 
cell. Hence, we will use the labels A (B) for even (odd) 
layers and the respective parameters. 

The full 0* theory for the magnetic fluctuations is 

= S2A + S2B + SiA + SiB + »5dyn + Sab- (7) 

The first two terms, S2A + S2B, contain the Gaussian 
description of the unfrustrated magnetism of the two 
tetragonal subsystems: 

S2A = 1 dT(fk\\ [(m^A + cifc|)(^^(fc||) 

n even 

+ riA^t'nA ■ (l>n+2,A , (8) 

the action S2B for odd layers is obtained hy B. Here 

0(fc||) is the order-parameter field after in-plane Fourier 

transformation, the momentum /c|| is now measured rel- 
ative to the ordering wavevector (tt, tt), and c is a ve- 
locity. Further, i]' represents the imfrustrated vertical 
second-neighbor coupling. From the microscopic model 
one reads off the bare value r]' J^z', however, as shown 
in Ref. 20, interaction effects cause a non-zero r/' even for 
vanishing bare Jzz, see Sec. VA. The mass is the 
control parameter of the zero-field transition; a negative 
(renormalized) mass leads to a magnetic condensate 

with {4>a) 7^ 0. (The coupling between (j)A and (pB will be 
discussed below.) The terms S^a, S^b contain the local 
quartic self-interaction. 



>4A 



E 

n even 



(9) 



Finally, 5dyn encodes the dynamics of the spin 
fluctuations:^ 

Sdyn = f dTd\ Y^^dJn ' iH X $nf ■ (10) 

In zero field, there is only a second-order time derivative, 
the dynamical exponent is 2; = 1, and the modes are 

triply degenerate. In contrast, in finite field we have z = 
2, and the modes are Zeeman-split according to w ^ w — 



ai?, where a = 0, — and the +, — modes correspond 
to (px i iipy (assuming the field to be in z direction). 

Let us now turn to the frustrated coupling between ad- 
jacent layers, described by Sab- To this end, we repeat 
the central symmetry argument for the perfectly frus- 
trated geometry, already given in Ref. 20. Apart from 
spin rotation and space inversion symmetry, the system 
is also invariant under QO-degree in-plane rotations, how- 
ever, the geometry dictates that this is accompanied by a 
relative sign change of the order parameter in two neigh- 
boring planes (!): 



-k^ , cj)n (-l)"ci„ 



(11) 



This symmetry strongly constrains the allowed inter- 
layer coupling terms in the case of perfect frustration. 
A general form of the inter-layer coupling, including 
quadratic and quartic terms, is 



Sab = j drd^kw ^ [n$n ■ $n+i + V^xky^ 



-1-1 



(12) 



(To avoid clutter of notation, we have assumed the cou- 
plings to be vertically unmodulated; a corresponding gen- 
eralization is straightforward.) Clearly, the k term ac- 
tually corresponding to unfrustrated inter-layer coupling 
- is incompatible with the symmetry (11), whereas the 
other terms are compatible. Thus, perfect frustration 
implies k = 0. The 77 term represents single-particle hop- 
ping in the presence of frustration, i.e., it vanishes for 
A;^; = or fcj^ = 0; the form of this hopping term can be 
directly obtained from expanding the tight-binding dis- 
persion on the bet lattice near in-plane wavevector (tt, tt). 
From the microscopic model (1) one reads off 77 oc J^. 

Further, wi and U2 represent density interactions be- 
tween adjacent layers. An important role in particu- 
lar in the ordered phase - is played by the ui term: in 
the presence of perfect frustration this is the leading cou- 
pling between adjacent layers at = which is sensitive 
to spin directions. Negative ui stabilizes collinear spin 
correlations between adjacent planes, while positive u\ 
favors orthogonal (f)A and (pB- Comparing with the mi- 
croscopic model (1), it is obvious that the ui term only 
arises at order J|, in fact ui cx —J'^/J^ see Sec. VA. 
(These statements hold in the absence of a bare Hcoii, 
otherwise ui oc JcoU dominates.) 



B. Full frustration: Z2 symmetry and Ising bond 
order parameter 

As the bilinear magnetic coupling between adjacent 

layers, • (pn+i, is suppressed by a prefactor of k^ky 
in the presence of perfect frustration, there is no linear 
coupling between the condensates {(I)a) and ((^^s) on the 
even and odd layers in an antiferromagnetically ordered 
phase. Instead, the dominant A-B coupling is given by 
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the biquadratic term oc ui, which wiU select colhnear or 
orthogonal correlations between the two condensates, but 
will always leave a Z2 degeneracy intact, corresponding to 
a spin inversion in every second plane. This Z2 symmetry 
corresponds to a true symmetry for the antiferromagnet 
on the bet lattice, and will be spontaneously broken in 
the ordered phase. 

For an undistorted lattice, it is then useful to introduce 
a local Ising order parameter which is conjugate 

to (pn ■ <Pn+i and lives at zero in-plane wavevector. ^' 
transforms as a singlet under SU(2) spin rotations, and 
hence can be described by an unfrustrated ^'^ theory. 
Assuming a single-dimer unit cell, a plausible form is: 

n 

+'7**Ti-l/2*n+l/2 + >S5r4 + (S^dyn, (13) 

where <S*4 is again a quartic self-interaction, and the dy- 
namic term »S*dyn contains a second-order time deriva- 
tive. The physical content of ^ is encoded in its interac- 
tion with (j), where the leading term is trilinear: 



(14) 



with A a coupling constant. (Additional couplings 
etc. do not modify the physics to be discussed below.) 

Eq. (14) shows that the condensation of induces an 
unfrustrated vertical hopping through the term A(\l/)(/>„ • 
(pn+i, i-e., it breaks the frustration. Microscopically, ^ 
condensation is equivalent to spontaneous bond order, 
modulating the vertical magnetic couplings within 
each unit cell, as illustrated in Fig. lb. Within a purely 
magnetic model, 4* can be understood as a singlet bound 
state of two </> quanta, implying that mip is essentially 
given by 2m^, plus a correction arising due to an attrac- 
tion or repulsion of (j) quanta from Sab (12)- 

In the presence of phononic degrees of freedom, bond 
order causes lattice distortions, hence ^ condensation is 
a structural phase transition. 4" may condense either 
uniformly or with a non-trivial modulation along the c 
axis. (In principle, the action could be dominated by 
phonon effects this will not be considered.) 

Let us now discuss the implications for the phase di- 
agram of the spin model on the ideal bet lattice: The 
full theory + + S^\s, admits two distinct scenar- 
ios: (i) A single transition driven by the condensation 
of (j) here, the coupling A generates a non-zero expec- 
tation value for \1/ as well, because the ui term in Sab 
leads to non-zero (</>„ • </)„+i). (ii) Two transitions: First, 
^ condenses, which modifies the quadratic part of the (j) 
action, relieving the frustration, and <p orders in a second, 
subsequent transition. For the microscopically relevant 
parameters, wc find that situation (i) is generically real- 
ized, see Sec. VC2. 



C. Canted magnetism near Hd 

Near the critical field Hd, only the lowest of the 
Zceman-split triplet modes is relevant for the low-energy 
behavior. This lowest mode is ^ = 4>x + i<Py, and the 
critical theory can be formulated using a single complex 
scalar field $ with canonical boson dynamics and a mass 
(i.e. chemical potential) ^ (iJci — -ff) • Condensation of $ 
breaks a U(l) symmetry and leads to spontaneous order 
perpendicular to the field direction: the transverse stag- 
gered magnetization is given by whereas the longitu- 
dinal uniform magnetization is ($*$). Thus, the physics 
near Hd is that of a dilute Bose gas.'' (Complications 
again arise from the additional Z2 symmetry, the detailed 
discussion will be given in Sec. VD.) 

The above symmetry analysis in terms of (j) continues 
to apply, with the change that 4> is now a two-component 
vector representing the transverse staggered magnetiza- 
tion. However, care is needed when associating effective 
with microscopic couplings, as e.g. 0^ now carries a uni- 
form magnetization and hence couples linearly to a field. 
(Technically, this arises because the degrees of freedom 
of the uniform magnetization are no longer gapped.) 



D. Canted magnetism near Hc2 

At high fields the ground state of the system is fully 
polarized, and the elementary excitations are bosonic 
spin-flip quasiparticles. Upon decreasing the field, those 
will condense at Hc2, leading to a canted phase which is 
continuously connected to the canted phase established 
above Hd- Hence, the order-parameter description is 
identical to the one near Hd- The order parameter is a 
canonical boson now with a mass (iJ — i?c2)- The 
transverse staggered magnetization is again given by ($), 
whereas the uniform magnetization is Mgat — 



E. Phase transitions 

Most of the quantum phase transitions discussed below 
are at or above their upper-critical dimension (which the 
exception of the Ising and 0(3) transitions in d = 2), thus 
the critical exponents are known. The shift exponent ip 
of the finite-temperature phase boundary is given by the 
product of correlation length and dynamical exponents, 
uz, if the QPT is below its upper critical dimension; oth- 
erwise it can be obtained from the temperature depen- 
dence of the Hartree diagram determining the mass shift 
of the order parameter,"' and is given by tp = z/{d+z—2) 
[z = 2 for a BEC transition). 

In this paper, we intend to estimate relevant energy 
scales and to obtain the overall behavior of observables, 
primarily at zero temperature away from the phase tran- 
sitions. To this end, wc employ bare as well as self- 
consistent perturbation theory. These methods may 



7 



break down near criticality, and we comment on this be- 
low. 

The finite-T regime close to the ordering temperature 

is more difficult: The interplay of the frustration-related 
order-from-disorder mechanism and the Mermin- Wagner 
theorem is delicate, see Sec. V E. A reliable treatment of 
the finite-temperature transitions, including an estimate 
of Tn, is beyond the scope of the paper. 



IV. BOND-OPERATOR THEORY 

For a quantitative study of the coupled-dimer Heisen- 
berg model (1) we apply the bond-operator approach of 
Sachdev and Bhatt,"^^ with extensions proposed by Ko- 
tov et al.^^ and by Sommer et alP'^ While this method 
can in principle be applied at finite temperatures,^^ we 
will restrict the explicit calculations to T = 0. 

In this section, we present the formalism for a situ- 
ation with equivalent dimers, i.e. a bet lattice with a 
single-site unit cell (with lattice sites denoted by i). The 
generalization to inequivalent dimers is straightforward, 
and we shall refrain from showing the lengthy equations. 



A. Paramagnetic phsise: Harmonic approximation 

The four states of a dimer i can be represented using 
bosonic "bond" operators {sl,^!^,} (a = x,y,z), which 
create the dimer states out of a fictitious vacuum. Ex- 



plicitly (and omitting the site index 



no), 



|a) = 4|0), where \s) = (|n) - |iT))/y2, \x} = 

(-im+iu))/A \y) =i(|tt)+iu))/Ai^) = (in)+ 

|it))/\/2. The Hilbert space dimension is conserved by 
imposing the constraint slsi+J2a Aa^ia = 1 on every site 
i. (In the presence of a Zeeman field, a rotated triplet 

basis with 4+ = (4 + it\y)/^, = - it\y)/^, 

Ag = Az useful as well.) 

The Heisenberg Hamiltonian (1) can now be formu- 
lated in terms of the bond operators {si,tia}, for details 
sec Rcfs. 28,32,34. To treat the paramagnetic phase, the 
following re-interpretation of the formalism is useful:^^ 
Starting from a background product state of singlets on 
all dimers, = Y\is\\0), the operators can be 

viewed as creating local triplet excitations in the sin- 
glet background. (Formally, this is achieved by setting 



1 which implies t 



tj^^.) The constraint 



then takes the form tla^ia < 1- (So far, the procedure 
is exact.) 

Upon expressing Ti. (1) in bond operators, products of 
two spin operators convert into terms with two, three, 
and four triplet operators, H = H2 +H4; the bi- 

quadratic spin term in TicoU (4) contains up to 8 triplets. 
For the ideal bet structure and an external field in z di- 



rection, the bilinear part W2 reads: 

qa 

(15) 

with a = +,0,-, a = -a, A^ = J + B^, = 2J'^^\ + 
'^Jizlqz and 

7g^l = (cosga;+cosgj^)/2, 

= cos(gx/2)cos(gj,/2)cosg'^. (16) 

The coupling constant appearing in A^^ B^ is J2Z = 
Jl^ + — Jl"^ — Jl^. The 3d momentum q runs over 
the Brillouin zone of the bet lattice, spanned by the 
primitive translations qi = (27r, 0, — tt), q2 = (0,27r, — tt), 
^3 = (0, 0, 2tt) in reciprocal space. 

The harmonic (or linearized) approximation consists 
in treating only the bilinear part Ti.2 of the Hamiltonian; 
both the hard-core constraint and the higher-order terms 
H3, Hi (given in App. A) are neglected. With the Bo- 
goliubov transformation t^^ = 
tonian TI2 can be diagonalizcd, with eigenvalues 



-VffT^ - „ the Hamil- 



and Bogoliubov coefficients 



aH 



Aq 



2uj, 



'qO 



Bq 



2u). 



qO 



(17) 



(18) 



In this linearized bond-operator theory, interactions be- 
tween the order-parameter fluctuations - represented by 
triplet quasiparticlcs are ignored. 

For \J2z\/J < 2 the dispersion minimum of the mag- 
netic modes is at q^^ = (tt, tt). There, ui^a is indepen- 
dent of qz due to frustration - hence effectively two- 
dimensional. Expanding near (tt, tt) yields: 



i^qO 



A+ ^ (^kj +r]kxky cos qz^ (19) 



where = q\[ — (t^iTt) and rj = —J2z/(2J'). The spin 
gap and velocity are given by A = J(J — 4 J') and 
c=^/jyj2. 

For larger unit cells, multiple triplon operators are in- 
troduced, and the Bogoliubov transformation needs to be 
performed numerically. 



B. Beyond the hcirmonic approximation 

Triplon interactions effects are important to lift the de- 
generacy of the dispersion along (tt, TT,qz). The most im- 
portant interaction correction arises from the hard-core 
constraint, which is conveniently implemented using an 
infinite on-site repulsion between the bosons:^^ 

Hu = uY^ 44/3*i«%, C/ ^ 00 . (20) 

iaj3 
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As proposed by Kotov et al."^^ , this hard-core term can 
be treated by a ladder summation of scattering diagrams, 
together with a self-consistent one-loop approximation 
for the self energy, Fig. 1 of Ref. 28. This method is 
also known as Brueckner approach in particle theory, 
the small parameter being the density of triplet bosons, 
which at zero temperature is given by ^ fS^ . 

Here, we have employed the formalism of Ref. 28 at 
finite fields (App. A). In the present problem, it is not 
sufficient to treat the quartic terms in 0.4 in a mean-field 
(Hartree-Fock) approximation: processes of second order 
in Ti.4 are needed to obtain the leading contribution cx 
to the unfrustrated second-neighbor hopping in z direc- 
tion. The calculations involve the self-consistent solution 
of integral equations and are performed numerically on 
lattices with up to 16^ sites. 

The approach of Ref. 28 has been shown to give 
results in good quantitative agreement with Quantum 
Monte Carlo and series-expansion methods, regarding 
e.g. phase boimdaries and magnetic excitations of the bi- 
layer Heisenberg model. Although the Brueckner method 
contains a re-summation of an infinite series of diagrams, 
it is not designed to capture critical behavior beyond 
mean- field. 



C. Ordered phases 

The bond-operator method can be generalized to mag- 
netically ordered states by taking into account the ap- 
propriate condensate. Technically, an expansion is then 
performed around a symmetry- broken product state IV'o) 
which replaces the singlet state. As demonstrated by 
Sommer et al.,'^^ a consistent description of the fluctu- 
ations is obtained by applying a harmonic approxima- 
tion after a rotation of the basis vectors in the four- 
dimensional Hilbert space of each dimer. 

Here we adopt the formalism of Ref. 29, briefly sum- 
marized in the following, to the bet lattice geometry. The 
rotated basis operators, replacing {sj,4a}' 



it 



ft 



1 





A2 


1 




VI + 


A2 


1 





_Ae^'5«^4 + 



(4 + i^Ay) 



it = it 



(21) 



where A and are condensate amplitudes. The role of 
the singlet product state is now taken by |?/'o) = Hi ^1 
For A = /z = we have the original "paramagnetic" bond 
operators, while A = 1, /i = describes a classical Necl 
state IV'o) and its local excitations (note that we have 
chosen the direction of the staggered magnetization to 
be in x direction). Finally, for /i = 1 and A ^ 00 the 
product state is the fully polarized state with all spins in 



z direction. Below, we will employ ordering wavevectors 
Q = (tt, tt, 0) or (tt, tt, tt), both describing states with fer- 
romagnetic correlations between 2nd vertical neighbors 
and reflecting the Z2 degeneracy w.r.t. the relative ori- 
entation of neighboring layers. 

The Hamiltonian can be re-written in terms of the 
{s,ta} operators; the corresponding lengthy expressions 
can be found in Ref. 29 and will not be reproduced here. 
The condensate parameters A and ji are determined by 
minimizing (^ol-H'IV'o); this can be shown to be equiva- 
lent to eliminating linear P terms in W. Then, as in the 
paramagnetic case, the ij^ can be treated as excitations 
on top of a background state, and s will be formally set to 
unity. Subsequently, the Hamiltonian admits a harmonic 
approximation, by only keeping the quadratic terms in 
ta ■ The resulting 7i2 is solved by a Bogoliubov transfor- 
mation, leading to new quasiparticles t^-q- For A = 1 and 
/i = 0, is equivalent to conventional linear spin- wave 
theory of a Neel-ordered antiferromagnet. 

In summary, the modified bond-operator approach in- 
terpolates between the triplon description of the param- 
agnet, spin waves of the antiferromagnet, and the flipped- 
spin quasiparticle physics of the field-polarized ferromag- 
net. Within the harmonic approximation, the phase tran- 
sitions turn out to be of second order, and the ordered 
phases have the correct number of Goldstone modes. We 
also note that an extension beyond the harmonic level is 
not obvious: Taking into account a hardcore repulsion 
of the ta as above leads to a violation of the Goldstone 
theorem. Therefore our quantitative calculations in the 
ordered phases below will be restricted to the harmonic 
level. This is expected to be a reasonable approximation 
away from the phase transitions, and is also qualitatively 
correct near phase transition above the upper-critical di- 
mension. However, we cannot capture critical behavior 
beyond mean-field, including possible logarithmic correc- 
tions occuring at phase transitions being at the upper- 
critical dimension (like a BEG transition in 2d). 



D. Observables 

In the harmonic approximation, static observables like 
magnetizations are calculated as expectation values with 
the ground state of H2 which is the vacuum of the 
Bogoliubov-transformed operators t^^. 

It may be tempting to calculate e.g. the T = uniform 
magnetization from the field dependence of the ground- 
state energy, M = —dEo/dH. However, this procedure 
is incorrect in magnetically ordered phases for the follow- 
ing reason: In an ordered phase, the condensate param- 
eters and hence the basis vectors are in general field- 
dependent, i.e., the expansion is done around a field- 
dependent product state. This implies that the quadratic 
part of the Hamiltonian will be field-dependent. In other 
words, in the harmonic approximation W2 corresponds to 
a different Hamiltonian for each field. Then, —dEo/dH 
contains, apart from the magnetization, an additional 



contribution arising from the field dependence of 7^2- We 
note that M = —dEo/dH is sometimes used in spin- wave 
theories for canted antiferromagnets - there, it yields in- 
correct results as well. However, the deviations from (Sz) 
are often small. 

Dynamic properties, like the cross section for inelas- 
tic neutron scattering, can be expressed in terms of the 
Green's functions of the triplon quasiparticles. Beyond 
the harmonic approximation, this route needs to be taken 
for static observables as well. 



E. Relation to single-boson description 

Let us close this section with comments on the relation 
between the advocated bond-boson approach and effec- 
tive theories based on a single-boson description often 
used in the field-driven cases. '^^ 

It is apparent, that near Hd only the t+ boson is low in 
energy (equivalent to the $ = <j)x+i4>y mode of Sec. Ill), 
hence is the only required degree of freedom in a crit- 
ical theory for the Hd transition. There are, however, 
a few subtleties: (i) Scattering processes between low- 
energy and high-energy bosons contribute to terms in the 
low-energy theory. This is in fact the case in the ideal 
bet lattice when J^z = in the Hamiltonian. Then, the 
bare boson dispersion is independent of Qz at wavevec- 
tors (7r,7r, and interactions are required to lift this 
degeneracy. As scattering processes between r+ bosons 
conserve particle number and hence are absent at T = 0, 
one may conclude that the dispersion remains flat. This 
is incorrect, because scattering processes with the high- 
energy particles r__o are no longer number-conserving, 
inducing a finite dispersion, (ii) To fully describe the con- 
densate with canted order throughout the phase diagram, 
a single complex parameter is insufficient, as clearly seen 
from the full bond-operator description (where both A 
and /I are needed). Using a condensate of t+ only is ap- 
propriate in the limit J' <C J, but this cannot describe 
states close to a quasiclassical Neel state. 

This concludes our description of the methods. 



V. CASE (A): EQUIVALENT LAYERS AND 
FULL FRUSTRATION 

This section discusses in detail the situation (A) with 
an ideal bet structure of dimers, while situation (B) - 
a structure with inequivalent layers and two dimers per 
unit cell - will be studied in Sec. VI. For the purpose 
of comparison with experiments, most numerical results 
shown below are for parameter values of possible rele- 
vance to the material BaCuSi206. 

The dimensional reduction observed in BaCuSi206, 
i.e., the 2d value of the shift exponent describing the 
BEG phase boundary, gc {H - Hd)"^ with V = 1, 
was originally discussed in terms of frustration only, i.e. 
scenario (A).^^ While the neutron scattering results of 
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FIG. 3: Diagrams occurring in the perturbation expansion 
of the order-parameter theory (7). The solid lines are (j> 
propagators (depending on in-plane momentum fcy and layer 
index n), the full circle is the local four-point vertex (oc uq), 
the cross is the frustrated inter-layer hopping (oc rjk^ky). a) 
Inter-layer density interaction [open circle, this includes both 
ui and U2 terms in Sab (12)], generated from a Uq process. 
For negativeui, the u\((l)n-<pn+i)'^ term leads to coUinear spin 
correlations in z direction, b) Unfrustrated second-neighbor 
vertical hopping (star, oc rj'), generated from interaction pro- 
cesses. This hopping is responsible for 3d behavior at lowest 
energies, irrespective of the inter-layer frustration, c) Addi- 
tional diagrams present in the ordered phase; the open square 
denotes the coupling to the condensate. The first diagram cor- 
responds to nearest-neighbor vertical hopping, the second is 
the leading "vertical" contribution to the free energy. 

Rcf. 21 point to a somewhat different origin of quasi-2d 
behavior, namely scenario (B) with inequivalent dimers, 
the physics of the ideal bet structure is a very interesting 
and instructive situation to study. Some of the results 
described below were reported by us in Ref. 20. 



A. Paramagnetic phase: Dispersions 

For one dimer per unit cell, there will be a single branch 
of magnetic triplet excitations in the zero-field singlet 
phase. Its in-plane dispersion caused by J' is unfrus- 
trated, but the c axis dispersion is heavily influenced by 
frustration: Away from the in-plane dispersion minimum 
at g|| = (tt, tt) it is dominated by the frustrated hopping 
between adjacent layers, which leads to a (bare) disper- 
sion of the form cos((7a;/2) cos(gj,/2) cos g^. This van- 
ishes at = (tt, tt), and only second- neighbor hopping 
contributes, with an additional dispersion Jzz cos2qz. 

A key question, also relevant for the critical behav- 
ior, is the fate of the vertical dispersion in a model with 
bare Jzz = 0. The arguments in Refs. 18 suggested 
that the dispersion would then be g^-independent along 
(tTjTt, g^), rendering the critical behavior at Hd truly 
two-dimensional. However, interaction processes invari- 
ably induce a symmetry-allowed effective vertical second- 
neighbor hopping^° of order Jz/ - this effect is the one 
responsible for the absence of true dimensional reduction 
(with the exception of the high-field situation H > Hc2, 
see Sec. VD3 below). 

Let us briefly repeat the perturbative arguments for 
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a) 



c) 



a. £)-"• (X 



FIG. 4: Diagrams in the Brueckner bond-operator approach. 
The sohd hnes are t^fa propagators, here depending on the fuU 
3d momentum q. a) Self-energy from the hard-core repulsion; 
the shaded square is the effective four-point vertex F (A5) 
obtained from a ladder summation, b) Leading self-energy 
from the quartic term TI4, with the circle corresponding to 
J4. c) Second-order self-energy in Ti.4 - its Til^ portion is 
required to obtain the leading term in the vertical dispersion 
on the ideal bet lattice, d) Second-order self-energy in Tisz, 
where the triangle is the three-point vertex of strength J^. 
The momentum space structure of the vertex suppresses this 
diagram at g]| = {■k,tt). 



the vertical dispersion. Within the effective order- 
parameter theory for <j), we need to look for processes 
which generate a contribution to the 77' term in Sa, 
Eq. (8), from the inter-layer interactions Sab, Eq. (12). 
The leading diagrams are in Fig. 3b, and are oc Ugry® 
and oc 2??^ . The microscopic identification rj (x Jz, 
1*1,2 oc suggests that both diagrams are of order . 
Higher-order processes will not change this result (except, 
perhaps, at a critical point, see Sec. VC below). How- 
ever, the order parameter theory misses interactions be- 
tween staggered and uniform magnetization fluctuations, 
as the latter are not contained in 5^. This becomes clear 
when discussing the same physics in the bond-operator 
language. Relevant interaction terms are the hard-core 
term Tiu (20) and the three- and four-point vertices aris- 
ing from J2, (A3), and (Al). (Note that the 
physics of Hsz and 'H4Z is absent from the effective order- 
parameter theory.) We will ignore effects of the in-plane 
part of Ti.4 beyond the Hartree contribution, because the 
dominant overall renormalization arises from Tijj and 
qualitative changes from 'Haw '^'^^ expected. Taking 
into account TLu yields a dispersion along (7r,7r,(7z) pro- 
portional to cos(2(7z). This is easily understood, as 
the self-consistent Hartree diagram. Fig. 4a, contains the 
process in Fig. 3b2. Now consider the 7i|^ portion of 
the diagram in Fig. 4c. Remarkably, it does generate 
a vertical dispersion proportional to (see App. A) ~ 
hence this the leading contribution to the interaction- 
generated vertical dispersion. We have therefore incorpo- 
rated this diagram into the Brueckner approach. Finally, 
there is Tizz- Fig. 4d shows the leading self-energy dia- 
gram, which, however, turns out to vanish at q\\ = (tt, tt). 
The contribution from T-L^z to the vertical dispersion is 
of order ; hence we will ignore Ti^z altogether. 

A full result for the dispersion is shown in Fig. 5. 
The parameter values are chosen with an eye towards 
BaCuSi206, i.e., we have tried to match the in-plane 
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FIG. 5: Triplon dispersion (solid) and boundary of the two- 
particle continuum (shaded), calculated using the Brueck- 
ner bond-operator approach (Sec. IV B and App. A) on an 
8^ lattice. The parameter values are J = 1, J' / J = 0.15, 
— J2z ~ Jaz ~ 0.2 J', Jsz = 0, H = 0. The inset shows an 
energy zoom into the almost flat vertical dispersion. For com- 
parison, the dotted line shows a dispersion calculated within 
the harmonic approximation, with parameter values chosen to 
match the Brueckner results at wavevectors (0, 0, 0), (tt, tt, 0), 
and (0,0, tt): J = 1.047, J'/J = 0.140, -J2z/J' = 0.196. 



dispersion averaged over the modes observed in Ref. 21, 
i.e., J'/J = 0.15. The vertical coupling is somewhat 
arbitrarily chosen as —J2Z = Jiz = 0.2J', giving a dis- 
persion at the unfrustrated point, i.e., along (0,0, gz), of 
(—0.096 Jcos ^2). The dispersion at the frustrated point 
is tiny, roughly (—2- 10^^" J cos 2^2), its bandwidth scales 
as as expected. 

In Fig. 5 we illustrate that the in-plane dispersion 
obtained from the Brueckner approach can be repro- 
duced using the harmonic bond-operator approximation 
describing non-interacting triplons with renormalized pa- 
rameters. (Similar observations were made earlier e.g. in 
the context of spin ladders. '^^) The leading renormaliza- 
tion is in J: triplon repulsion pushes the dispersion to 
higher energies, hence the "harmonic" J is larger than 
the "true" (Brueckner) J. Of course, the renormaliza- 
tion in general depends on temperature, magnetic field 
and other parameters. Using the zero-field renormalized 
parameters gives a reasonable account of the field depen- 
dence of the spin gap, and overestimates Hd only by a 
few percent (compared to the Brueckner approach). In 
the ordered phases and in the more complicated case of 
inequivalent layers (Sec. VI), where a Brueckner calcula- 
tion is no longer feasible, we will exploit this fact, i.e., we 
will work with the harmonic approximation and param- 
eter values chosen to match experimental data. 

Moving closer to the zero-field transition by increasing 
J'/J, the overall triplon bandwidth increases, see Fig. 6. 
The interaction-generated bandwidth Ez along (tt, tt, Qz) 
strongly increases, as the relevant energy denominator in 
the fourth-order expression is given by the third power of 
an averaged triplon energy. The induced bandwidth will 
also depend on temperature, but this effect is exponen- 
tially suppressed for T < A. From the numerical results, 
we are not able to track Ez close to the phase transition, 
due to discretization errors on the finite lattice. Finally, 
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upon applying a field, the triplon modes simply split as 
expected, Fig. 7. 

After having established the vertical second-neighbor 
coupling, it is worth discussing the first-neighhoT corre- 
lations. In the field-theory language, these arise from the 
Ml term in Sab (12) and are of orthogonal (coUinear) 
type for ui > (ui < 0), as the ui term is insensitive to 
the sign of 4>a,b- These correlations can be seen in four- 
point spin correlators only. The peculiar momentum- 
space structure of the i] term implies that two-point 
spin correlators between even and odd planes are strictly 
zero within the order-parameter theory in its paramag- 
netic phase. (The microscopic bond-operator approach 
shows weak ferromagnetic correlations between neighbor- 
ing planes, which do not become long-ranged near the 
ordering transition.) 

Microscopically, the ui term may originate from a bi- 
quadratic spin coupling as in Ticoii (4). As above, one has 
to ask about the value of ui in the absence of a bare Jcoii- 
In fact, ui^2 will be interaction-generated as well, with the 
leading diagram cx shown in Fig. 3a. Care needs to 
be taken when contracting two uq vertices because of the 
vector structure of the interactions. A straightforward 
analysis shows ui, U2 < with M2/U1 — (A^ + 4)/4, where 
N is the number of order parameter components, ui < 
means that collinear correlations are favored - this is con- 
sistent with the quasiclassical order-from-disorder mech- 
anism operative inside the ordered phases, see Sec. VB. 
(Note that U2 acquires additional contributions cx Jz near 
Hci,2 from a mean-field-like interaction of uniform mag- 
netizations.) 

B. Antiferromagnetic phase: Order from disorder 

We now turn to the physics inside the antiferromagnet- 
ically ordered phases; the discussion of the phase transi- 
tions will be postponed to the following subsections. 

As mentioned in the introduction, for classical mo- 
ments on the bet lattice, antiferromagnetic layers are 
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FIG. 6: Triplon dispersion (solid) and two-triplon continuum 
(shaded) as in Fig. 5, but with parameter values closer to the 
zero-field transition: J = 1, J'/ J = 0.25, — J2z — Jaz = 0.2J', 
Jsz = 0, -H" = 0. Compared to Fig. 5, the bandwidth along 
(tt, 7r,gz) is significantly enhanced. 
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FIG. 7: Triplon dispersion as in Fig. 5, but now in a finite 
field; J = 1, J'/ J = 0.15, -J2z = Jiz = 0.2 J', Jsz = 0, 
H = 0.4. 

completely decoupled (provided J^z — 0). 



1. Spin-wave theory 

The effect of quantum fluctuations on the frustrated 
bet antiferromagnet has been initially studied using semi- 
classical spin- wave theory, '^^"^^ and we briefly summarize 
the results here. 

The simplest calculation, for a model with spins S 
on sites of a bet lattice, assumes a spiral order with 
wavevcctor (7r,7r,(5z). In linear spin-wave approxima- 
tion, the spin-wave spectrum is degenerate along the line 
(7r,7r, g^), signaling frustration. However, the ground 
state energy contains an inter-layer contribution of the 
form {— J^S/ J' cos^ Qz), thus favoring collinear order 
with Qz = or Qz = tt - this term arises from the zero- 
point energy of high-energy spin-wave modes. Including 
spin-wave interactions to order 1/S removes the degen- 
eracy in the dispersion along (7r,7r,(7z), consistent with 
the assumed ordering wavevector - the dispersion is pro- 
portional to ^/S\Jz sm{qz/2)\, i.e., corresponds to vertical 
flrst-neighbor hopping."*" Thus, this calculation predicts 
ferromagnetic order between second-neighbor planes, and 
collinear order between neighboring planes. A subse- 
quent, more detailed, calculation*^ relaxed the assump- 
tion of a single ordering wavevector and analyzed more 
general ordered states to high orders in 1/S. As a re- 
sult, for S — 1/2 a state with antiferromagnetic order 
between second-neighbor planes and collinear order be- 
tween neighboring planes was found be selected by energy 
contributions of order J^/{JSf and J^/{J^S). Thus, 
fully 3d order is stabilized within the stacks of even and 
odd planes. The residual Z2 degeneracy specifying the 
relative orientation of the "even" and "odd" order pa- 
rameter (Sec. IIIB) is left intact and hence spontaneously 
broken in the ordered state. 

We note that the tendency towards collinear correla- 
tions can be mimicked^* by a biquadratic term in the 
Hamiltonian of the form (4), with Jcoii c>c — J|/J' - we 
will exploit this in the triplon bound-state calculation. 

One may also consider order-from-disorder on the 
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bet lattice in the different situation with non-vanishing 
second- neighbor vertical coupling J^^. Then, 3d order 
in each of the tetragonal subsystems is already estab- 
lished at the classical level, but the two order parameters 
are decoupled. Quantum fluctuations are only needed 
to produce the colHnear coupling (oc Jf/J'), reducing 
the degeneracy to the residual Z2 symmetry. (This lat- 
ter order-from-disorder mechanism is similar to the one 
in the much-studied Ji — J2 Heisenberg model on the 
square lattice, in the limit of large J2.) 

2. Field theory 

It is instructive to re-phrase the order-from-disorder 
mechanism in terms of the order-parameter field theory 
of Sec. Ill, with the relevant diagrams shown in Fig. 3. 

First, the magnetization orientations between adjacent 
layers are determined by the ■(f'n+i)'^ term. Indeed, 

we have argued above for ui < which favors coUinear 
correlations. Second, the type of 3d ordering of second- 
neighbor layers is determined by the Tj^cp^ • 

term. 

Third, the coupling to the condensate allows an effec- 
tive nearest-neighbor hopping in z direction, which was 
symmetry-forbidden in the paramagnetic phase. The di- 
agrams. Fig. 3c, show that the order-from-disorder con- 
tributions to the vertical dispersion and to the free energy 
scale with the square and the fourth power of the order 
parameter, respectively. 

3. Bond operators 

In the Brueckner bond-operator approach, we can fol- 
low the properties in the paramagnetic phase up to the 
phase transition. From the single-particle dispersion. 
Fig. 5, we can read off t]'; we find rj' < 0, i.e., a dis- 
persion minimum at (7r,7r, 0) and (tt, 7r,7r). Triplon con- 
densation at this wavevector implies ferromagnetic orien- 
tation between second-neighbor layers. This may appear 
inconsistent with the above-mentioned result of Ref. 41, 
but the latter strictly only applies in the semiclassical 
limit. The coUinear inter- layer correlations are contained 
in four-point triplon correlators. In those cases where we 
find a triplon bound state (Sec. VC2 below), its inter- 
nal structure implies coUinear correlations and its disper- 
sion has a minimmn at wavevector (0,0,0). Condensa- 
tion of this bound state implies uniform bond order (see 
Sec. IIIB) and again ferromagnetic orientation between 
second- neighbor layers. [The ordering pattern of Ref. 41 
would correspond to condensation of ^' at (0, 0, tt).] 

Inside the ordered phases, our bond-operator calcula- 
tions are restricted to the harmonic approximation. Here, 
the dispersion degeneracy along (tt, tt, g^) is not lifted: 
there are too many Goldstone modes. Nevertheless, we 
expect energy-integrated properties like magnetizations 
to be semi-quantitativcly correct. Some results in the 
field-induced canted phase are shown in Fig. 9 below. 



C. Zero- field quantum phcise transition 

The quantum paramagnet can be driven into an anti- 
ferromagnetically ordered state at zero field by increas- 
ing J'/ J. Experimentally, the application of pressure 
changes J and J' via a modification of bond lengths and 
angles. The coupled-dimer material TlCuCls is driven 
into an ordered state upon application of pressure, hence 
J'/ J is increased here; for BaCuSi206 systematic pres- 
sure studies have not been performed to our knowledge. 

In an unfrustrated system, the zero-field ordering tran- 
sition breaks the SU(2) spin symmetry and has dynamic 
exponent z = 1. In a spatially anisotropic, i.e., layered, 
system with a 3d scale E^, the quantum critical behavior 
is 2d for energies or temperatures above while it is 
3d below E;^. Importantly, E^ is essentially given by the 
bandwidth of the vertical dispersion near the ordering 
wavevector. 

For the frustrated bet system, several complications 
arise which we briefiy address in the following: There is 
the additional Z2 symmetry to be broken, and the ver- 
tical dispersion consists of a bare frustrated part and an 
interaction-generated unfrustrated part. 

1. Interaction-generated dispersion 

If the bare Jzz vanishes, then a vertical dispersion 
along {ir.TT.qz) (determining the 3d crossover sCf lie) only 
arises from interactions. The leading term is of second 
order in a four-point vertex (Sec. VA), it strongly de- 
pends on the distance to criticality because the spin gap 
enters the denominator of the perturbative expression. 

The simplest perturbative estimate is reliable in the 
paramagnetic phase for not too large J^, but may fail at 
criticality due to a non-trivial flow of the four-point ver- 
tex. This happens at or below the upper critical dimen- 
sion d+ and is signaled by infrared divergencies in the 
next-order diagrams. The renormalization group thus 
has a rather interesting structure (assuming a vanish- 
ing bare rj'): one-loop accuracy is sufficient to study the 
flow of the four-point vertex, but the strongly relevant 
two-point vertex rj' is only generated at two loops. (A 
similar situation has recently been studied in the context 
of spin chains with a frustrated inter-chain coupling.*^) 
Here we shall not perform a consistent two-loop renor- 
malization group treatment; we expect that the scaling 
Ez oc continues to hold for small Jz, albeit with non- 
trivial finite-temperature corrections to the prefactor in 
the quantum critical regime. (The scale is expected 
to evolve continuously into the scale of the ordered 
phase. Sec. VB, which decides between ferro- and anti- 
ferromagnetic ordering between second- neighbor layers.) 

We also note that the bond-operator calculation par- 
tially accounts for the above renormalization effects due 
to the self-consistent resummation of vertex diagrams. 
(For the field-driven case, such a strong renormalization 
of 77' does not arise, as the interaction processes nec- 
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essarily involve the high-energy triplon branches, due 
to particle-number conservation within the low-energy 
branch.) 

2. Bound state condensation and split transition 

As discussed in Sec. IIIB, the Z2 symmetry can be 
broken via condensation of the bond order parameter ^' 
before magnetic order occurs. 5" corresponds to a sin- 
glet bound state of two triplons."^^ While the transition 
to bond order may be studied using a mean-field ap- 
proach, the bound-state dynamics requires more effort. 
As in Ref. 37, we have solved a Bethe-Salpeter equa- 
tion for two-triplon bound states. The required input 
encompasses the single-triplon dispersion - to be taken 
from the Brueckner bond-operator calculation - and the 
triplon four-point scattering vertex. As self-consistency 
in the two-particle sector is beyond reach, it is reasonable 
to work with a bare scattering vertex, with contributions 
from 7^4 and Tiu- However, from our analysis it is clear 
that the effective biquadratic inter-layer interaction, rep- 
resented by ui{(f>n ■ (f>n+i)'^ in the field theory, is relevant 
for triplon attraction. We therefore include by hand TicoU 
(4) which mimicks the effect of the coUinear interaction; 
details are relegated to App. B. 

Interestingly, for triplon attraction the couplings 
(— Jcou) and {+J4z) act in a very similar way: both can 
cause inter-layer binding. While the triplon binding ef- 
fect of Ticoii is easily understood in the order-parameter 
language, this is more subtle with 7^42. First, we ob- 
serve that the relevant part of 7i4z involves the uni- 
form (as opposed to staggered) degrees of freedom of 
the dimers. Second, ^' condensation represents bond or- 
der arising from quantum-mechanical singlet formation, 
which is common for frustrated spin-1/2 systems.'^* (Note 
that the attractive force is linear in J^z-) 

In zero field, we find a singlet bound state below the 
two-particle continuum only for sufficiently strong attrac- 
tion (J42 — Jcou)- This is plausible, as small attraction 
in 2d causes only an exponentially shallow bound state, 
while in 3d a finite attraction is needed for binding - 
this simply follows from the properties of one particle 
with quadratic dispersion moving in a potential well. The 
wavefunction of the lowest bound state changes sign un- 
der 90-degree in-plane rotations of the internal coordi- 
nate; from the symmetry considerations in Sec. Ill we 
conclude that this bound state corresponds to a quan- 
tum as anticipated. The bound-state dispersion has its 
minimum at (0,0,0) and a bandwidth along {0,0, qz) of 
J^, inherited from the two-triplon continuum. A conden- 
sation of this singlet bound state corresponds to the Z2 
symmetry breaking advocated in Sec. IIIB. The bound 
state is unaffected by a Zeeman field, hence for fields 
larger than the triplon binding energy the bound state 
ceases to exist below the two-particle gap. 

In conclusion, both of the following scenarios may be 
realized (Fig. 8): (i) Two transitions: Upon increasing 
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FIG. 8: Schematic T = phase diagrams for the bet lattice 
coupled-dimer model, for the situation of equivalent layers 
with perfect frustration, a) Pressure tuning - here generically 
two scenarios are possible (where the second one applies for 
the parameters relevant to BaCuSi2 06). b) Field tuning. 



J'/J, there is first a quantum Ising transition with z = 1 
where bond order is established through condensation of 
^f. This transition is asymptotically 3d: the vertical "if 
dispersion is of order J* along {0,0, Qz). Subsequently, 
magnetic order will be established in a second transition 
where (f) condenses and the SU(2) symmetry is broken. 
The second transition has z = 1 and is of conventional 
0(3) type, in the sense that frustration is removed by 
the 5* condensate, i.e., by the term A(^')0„ • [see 
Eq. (14)]. The distance of the transitions is determined 
by the triplon binding energy, (ii) One transition: the 
condensation of (j) breaks the SU(2)xZ2 symmetry in one 
step. This transition is again asymptotically 3d; a de- 
tailed study of the critical behavior shall not be under- 
taken here. 

For the parameter values of possible relevance to 
BaCuSi206, we see no indication for bound states; hence 
case (ii) applies. Unfortunately, the numerical Bethe- 
Salpeter calculation suffers from severe finite-size effects, 
therefore we cannot reliably resolve small binding ener- 
gies. We also note that the inclusion of lattice effects 
can modify the behavior: for instance, bond order could 
occur far before magnetic order sets in - this applies not 
only to r = but also to the finite-T transitions. 



3. Anisotropies and further perturbations 

Perturbations beyond the model Hamiltonian (1) can 
modify the asymptotic critical behavior: Those include 
magnetic anisotropies of Dshyaloshinski-Moriya or dipo- 
lar type, as well as coupling to nuclear spins. In partic- 
ular, anisotropy terms are relevant perturbations to the 
0(3) critical points, rendering the asymptotic critical be- 
havior Ising-like. 
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FIG. 9: a) Triplon energy gaps and b) uniform (solid) as well 
as staggered (dashed) magnetization, calculated as function 
of the external field H using bond operators in the harmonic 
approximation. Parameter values describe the ideal bet lat- 
tice: J = 4.66 meV, J' = 0.5 meV, -J2z = 0.1 meV. 



D. Field-driven quantum phase transitions 

1. Lower critical field Hd 

The ground state of the quantum paramagnet is unaf- 
fected by applying a Zeeman field, however, the triplon 
excitations will split, leading to a spin gap decreasing 
with applied field. At at critical field Hd, given by the 
zero-field gap, the lowest triplon mode becomes soft, and 
antiferromagnetic order perpendicular to the field is es- 
tablished. The phase transition is in the dilute Bose gas 
universality class, with dynamical exponent z = 2, and 
breaks a U(l) symmetry.^ 

As above, the energy scale Ez for the dimensional 
crossover between 2d and 3d is given by the vertical band- 
width of the single-particle dispersion near the minimum 
wavevector. In bet lattice coupled-dimer model with per- 
fect inter-layer frustration, this scale is again given by 
the sum of the bare vertical second-neighbor hopping, 
Jzz, and the corresponding interaction- induced contribu- 
tion - for the latter the perturbative estimate J^/J^ is 
now reliable even at the transition, because the relevant 
T = interaction processes involve high-energy triplons. 
At the transition, U(l) and Ising symmetries are broken 
simultaneously. Fig. 8b, as singlet bound states are ir- 
relevant in finite field. The phase transition is still of 
the dilute Bose gas type, i.e. with z = 2 and mean-field 
exponents, supplemented by logarithmic corrections in 
the 2d regime. (Magnetic anisotropics beyond the model 



(1) can modify the critical behavior at lowest energies as 
noted above; a Dshyaloshinski-Moriya interaction may 
even smear out the phase transition.) 

For numerical calculations we resort to the harmonic 
bond-operator method. To this end, we employ param- 
eters J = 4.66 meV, J' = 0.5 meV, -Ja^ =0.1 meV, 
which roughly reproduce (within the harmonic approxi- 
mation) the BaCuSi2 06 mode dispersions (averaged over 
the multiple modes). Results for mode gaps and magne- 
tizations are shown in Fig. 9 (here Q = (tt, tt, 0)). 

2. Phase boundary: 3d critical behavior at elevated 
temperature ? 

While inter-layer frustration leads to well-defined 2d 
quantum critical behavior at intermediate energies or 
temperatures, one can ask whether 3d critical behavior 
(with associated exponents) is restored at higher energies 
- this was proposed on the basis of the experiments on 
BaCuSiaOg in Ref. 16. 

While there is a sizeable vertical dispersion of the mag- 
netic modes away from in-plane wavevector (tt, tt), which 
may be associated with 3d behavior, the thermodynamics 
is (to leading order) determined by the density of states 
of all magnetic modes. Here the alternating sign of the 
inter-layer hopping (as function of in-plane momentum) 
turns out to be crucial, because the 3d-like contributions 
to the density of states tend to cancel. As a result, the 
density of states does never cross over to the power law 
characteristic for 3d - this applies to both zero-field and 
finite-field cases. 

Turning to the location of the boundary of the ordered 
phase: this can be estimated within the order-parameter 
field theory by calculating the finite-temperature mass 
correction, i.e., the temperature dependence of the order 
parameter self-energy, on the paramagnetic side. For a 
transition above the upper critical dimension, the lowest- 
order estimate given by the Hartree diagram is usually 
sufficient; what enters the Hartree diagram is precisely 
the density of states of the magnetic modes. Performing 
the integral numerically nicely shows a crossover from 
T.{T) - E(0) cx r3/2 for T below the tiny to E(T) - 
E(0) oc T above Ez, but no further crossover to another 
well-defined power law. 

This strongly suggests that experimentally observed 
deviations from the 2d critical power laws (above 1 K) 
have nothing to do with 3d critical behavior, but instead 
indicate that one leaves the critical regime (in the sense 
that the correlation length is no longer large, or that the 
density of triplons in the ground state is no longer dilute) . 

3. Upper critical field Hc2 

For large external fields, H > Hc2, the ground state of 
the coupled-dimer model is fully polarized, and a quan- 
tum phase transition to a canted state occurs at Hc2 
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which is in the dilute Bose gas universahty class as well. 

In contrast to the paramagnetic low-field phase, the 
wavefunctions for the ground state and the one-particle 
excitations above Hc2 are exactly known: A Bloch wave 
of one spin flip on top of the ferromagnetic background 
is an exact eigcnstatc of Ti. This has interesting conse- 
quences: although these spin-flip particles have a hard- 
core interaction, scattering processes are rare at low T, 
as the equilibrium particle density vanishes as T ^ 0. 
Thus, the interaction corrections to the bare dispersion 
are exponentially suppressed above i?c2-^^ 

Hence, in a model without bare Jzz, the vertical mode 
dispersion vanishes identically at T = for H > , and 
will only be induced by thermal processes in the quantum 
critical regime. This suggests that the asymptotic critical 
behavior is two-dimensional.^^ 



E. Classical phase transition and Neel temperature 

Given the order-from-disorder mechanism, which gen- 
erates an effective 3d coupling proportional to the square 
of the order parameter itself, and given the Mermin- 
Wagner theorem, the obvious question about the nature 
of the finite-temperature transition arises. The issue ap- 
pears particularly relevant for 0(3) symmetry, i.e., in the 
zero-field case: Starting in the ordered phase at low T, 
the system has a robust (effective) vertical coupling of or- 
der Jz, and one would hence predict a Neel temperature 
Tn roughly given by J' / \n{J' / J^) (in the limit of small 
Jz, and J' is the relevant in-plane coupling constant). '^^ 
However, in the high-temperature paramagnetic phase 
the vertical coupling is tiny, i.e., only given by the ef- 
fective Jzz oc Jz- Upon cooling from high T this results 
in an instability of the paramagnetic phase at a tem- 
perature which is smaller by a factor of 4 compared to 
the above estimate. This argument indicates that the 
finite-temperature transition could be discontinuous, but 
precursor effects of the order-from-disorder mechanism 
may counteract, rendering the transition continuous. A 
reliable self-consistent treatment of this problem appears 
difficult and is beyond the scope of the paper. 

In any case, it is clear that the Neel temperature will 
strongly depend on the frustrated inter-layer coupling Jz ■ 
This applies not only to the situation with 0(3) symme- 
try, but also to the 0(2) case: here, the ordering temper- 
ature without vertical coupling does not vanish, but is 
given by the Kosterlitz-Thouless temperature Tkt, and 
Jz leads to a logarithmic enhancement of Tn compared to 
Tkt- If the inter-layer coupling is the sum of an unfrus- 
trated and a frustrated part, then both will contribute to 
Tn - this is likely the scenario relevant to BaCuSi206. 

VI. CASE (B): INEQUIVALENT LAYERS 

We now turn to the case where the ideal bet struc- 
ture is distorted such that the unit cell is enlarged to 



contain multiple dimers. This modification is suggested 
by the results of zero-field inelastic neutron scattering 
on BaCuSi206, where indications for at least two modes 
were found. 

As discussed in Sec. II B, we restrict om* calculations to 
a situation with two types of inequivalcnt layers, A and 
B, which are stacked in an alternating fashion. (More 
complicated symmetry-breaking patterns only lead to 
quantitative modifications.) We will consider both the 
cases of perfect and imperfect frustration, the latter mod- 
eled by unequal vertical couplings along the two diago- 
nals, see Sec. II B. 

The different microscopic couplings Ja,b and g in 
the even and odd layers then translate into unequal order 
parameter masses m^A,B and velocities ca.b in Eq- (8), 
whereas imperfect frustration implies finite k in Eq. (12). 
(In the following we will assume m^A < Tn^ b ■ ) 

Numerical results will be obtained from the bond- 
operator theory in the harmonic approximation, applied 
to both the disordered and ordered phases. For sim- 
plicity, we will not take into account any effects beyond 
quadratic terms; this implies that we miss, e.g., effects 
of the coUincar coupling between the two condensates. 
This will affect some low-energy properties and will be 
noted below, but the gross features of the results can be 
expected to be correct. 



A. Paramagnetic phase: Dispersions and reduced 
dimensionality 

In the situation of two positive triplon masses, the ther- 
modynamic and magnetic properties of the system will 
be governed by the smaller of the two gaps. The primary 
signature of the doubled unit cell is the presence of two 
triplet modes at fixed wavcvector (instead of one). Let 
us focus on the dispersion of these modes, in a regime of 
small to moderate vertical couplings Jz (neglecting the 
tiny Jzz)- Then, the in-plane dispersion is conventional 
and dominated by the J, J' values. At fixed q\\. the en- 
ergy of the two modes differs by Aw,. For Jz <C Awg, 
the vertical dispersion arises only in second order in Jz 
(!), i.e., cx J| cos(2(j2)/Acjg. In contrast, for Acjg 
(which includes the case of equivalent planes) the vertical 
dispersion is oc Jz as usual. 

Hence, inequivalcnt layers present an alternative mech- 
anism for dimensional reduction, which moreover is ef- 
fective for all in-plane wavevectors and in the absence of 
perfect frustration. Explicitly: for imperfect frustration, 
the relevant 3d energy scale [being the bandwidth along 
(7r,7r, g^)] is now oc J| for small Jz- Note that perfect 
frustration results in a prefactor of cos^(gx/2) cos^{qy/2) 
in front of the bare vertical hopping, which leaves us with 
the interaction-induced contribution to the vertical cou- 
pling ai Qx = TT or Qy = TT, which is still of order J^. 

Let us make these statements quantitative, by fitting 
the BaCuSi206 neutron scattering data of Ruegg et al."^^ - 
As we model only two inequivalcnt planes, we take the 
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FIG. 10: Zero-field triplon dispersion for inequivalent layers 
with perfect frustration, calculated using the harmonic bond- 
operator approximation. Parameter values are: Ja ~ 4.27 
meV, Jb = 5.04 meV, J'a ^ J'g = 0.5 meV, -J2z = = 
0.1 meV. Compared to the equivalent- layer case, the vertical 
dispersion at the unfrustrated point, i.e., along (0,0, (?z), is 
both quenched and modified in shape, from oc Jz cos to 
oc ±jf cos 2^2. 
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FIG. 12: Bandwidths of the vertical dispersion along (tt, n, q^) 
and (0, 0, qz) as function of the average Jz, for equivalent lay- 
ers (J = 4.66 meV) and perfect frustration (solid, squares), in- 
equivalent layers {Ja = 4.27 meV, Jb = 5.04 meV) with per- 
fect frustration (dotted, triangles), imperfect frustration with 
Jzi/Jz2 = 3 (dashed, circles), and imperfect frustration plus 
vertical modulation with JzAi/JzA2 = 11/9, JzBi/JzB2 ~ 19, 
JzA ~ JzB (dash-dot, no symbols, see also Fig. 18) below. In 
all cases, J' = 0.5 meV. 



neutron peaks with smallest and largest energies and find 
the (effective) coupHngs Ja = 4.27 meV, Jb = 5.04 meV, 
J'a = J'b — 0-5 nicV, being essentially identical to the 
ones estimated in Ref. 21. Not much is known about the 
values of the vertical couplings in BaCuSi206. The sim- 
plest assumption is perfect frustration - a sample result 
for the zero-field triplon dispersion for J^ = 0.1 meV is 
shown in Fig. 10. The vertical bandwidth along (0, 0, Qz) 
is roughly 0.04 meV: compared to the equivalent-layer 
case, it is suppressed here by a factor of 5, and it fol- 
lows cos2(72 instead of cosg^ (!). The bandwidth along 
(7r,7r, g^) vanishes due to the perfect frustration. (The 
tiny contribution is not captured in the harmonic 
approximation.) For comparison, in Fig. 11 we also 
show dispersions in a situation with imperfect frustra- 
tion, with couplings along the two inequivalent diagonals 
of strength J^i = 0.15 meV, Jz2 — 0.05 meV. Here, the 
bandwidth at the unfrustrated point is as before, but 
along {TT,Tr,qz) we now observe a dispersion of width 
0.025 meV. 



In Figs. 12,13 we show vertical bandwidths as func- 
tion of Jz for the above values of in-plane coupling. 
For Jz <C Acjg, these bandwidths scale as {JzAi — 
JzA2){JzBi ^ JzB2) / ^<-^q at the frustrated point and 
{J zAi-\- J zA2){J zBi+ J zB2) I ^i^q at the unfrustratcd point 
[with the conventions of Eq. (6)]. In Fig. 13 the vertical 
bandwidths are shown as function of the two-dimensional 
in-plane momentum, for the parameter values of Figs. 10 
and 18 (below, this case features inequivalent vertical di- 
agonals). 

B. Antiferromagnetic phases: Frustrated proximity 

effect 

The antiferromagnetic phases in the inequivalent-layer 
case, both at zero and finite field, can be nicely discussed 
using the order-parameter field theory. Due to the un- 
equal masses for the A and B subsystems, a situation 
with m^A < 0, m^B > can occur. Then, order in the 
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FIG. 11: As in Fig. 10, but now for inequivalent layers with 
imperfect frustration. Parameter values are: J a = 4.27 meV, 
Jb = 5.04 meV, J'a = J'b ^ 0.5 meV, Jzi = 0.05 meV, 
Jz2 = 0.15 meV - the latter two values are the couplings 
along the two inequivalent Jz diagonals, see Fig. lb. 
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FIG. 13: Vertical dispersion of the lowest mode, 
uj{qx, qy,Ti /2) — uo{qx, qy, 0), as a function of in-plane momen- 
tum for inequivalent-layer situations, a) Perfect frustra- 
tion, parameter values as in Fig. 10. b) Imperfect frustration 
plus vertical modulation, parameter values as in Fig. 18) be- 
low. 
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FIG. 14: Schematic T = phcise diagrams for the coupled- 
dimer model with inequivalent layers and perfect frustration. 
Due to the absence of a proximity effect, there is no coupling 
between the condensates on the subsystems A and B (on even 
and odd layers), leading to two magnetic transitions, a) Pres- 
sure tuning, b) Field tuning. 

A subsystem is established. 

Usually, one would expect that A order induces order 
in the B subsystem as well, due to a proximity effect. 
Technically, a proximity effect arises from a linear cou- 
pling between (j^A and 4>b [the k term in Sci,ab (12)]. 
However, as detailed in Sec. Ill, this term is forbidden in 
the presence of perfect frustration, and hence the prox- 
imity effect is absent due to frustration! The only effect 
of the A condensate on the (t)B fluctuations is via the 
quartic terms which cause an anisotropy, see below. 

The absence of a proximity effect also implies that 
the (strong) order-from-disorder mechanism discussed in 
Sec. VB is only present in phases with order in both A 
and B subsystems. 



is larger than binding energy of two A and B triplons 
(which is set by J42), then the first transition is the mag- 
netic ordering of (j)A- This transition is in the standard 
0(3) universality class and has z = 1. The dimensional 
crossover between 3d and 2d critical behavior is set by 
ry^ ~ J^, as discussed in Sec. VC. 

Ordering of A introduces a mass anisotropy into the 
action for (j)B through the term ui{{4>a) • 0b)^, which 
is of easy-plane (easy-axis) type for ui > (ui < 0). 
Then, further increasing J'/J will lead to ordering of 
<f>B, with a quantum phase transition in the 0(2) (Ising) 
universality class with z ~ \. Recall that wi < is the 
generic situation in our coupled-dimer model, leading to 
a secondary Ising transition. The crossovers near the 
second transition are somewhat more complicated: with 
increasing energy, we go from 3d Ising to 2d Ising to 2d 
Heisenberg, with the two crossover scales set by ry^ ex 

and ui {4>a)'^ « J|. 

For imperfect frustration there is only one magnetic 
transition in the 0(3) universality class, and the sec- 
ondary transition is smeared out. The energy scale rele- 
vant to the dimensional crossover from 3d to 2d is now de- 
termined by the "bare" triplon dispersion near the mini- 
mum wavevector (tt, tt, 0) (i.e. interaction effects are sub- 
leading). As shown in Sec. VI A, this scale is proportional 
to (i.e. much larger than in the perfectly frustrated 
case). 

E. Field-driven phase transitions 



C. Phase diagrams 

Given the absence of a proximity effect in the perfectly 
frustrated case, two distinct phase transitions will occur 
for the A and B subsystems. Depending on the rela- 
tion between the parameters Ja,b, J a b^ different phase 
diagrams may be reaHzed, see Figs. 14,15. 

There will be primary transitions between a disordered 
and a partially symmetry-broken phase, and secondary 
transitions between partially and fully symmetry-broken 
phases. Both types transitions extend to finite temper- 
atures and will be visible as singularities in thermody- 
namic quantities. Breaking the frustration will turn the 
secondary transitions into crossovers, see Sec. VIE for 
explicit numerical results. 

We believe that the presence or absence of these sec- 
ondary transitions allows a clear-cut experimental dis- 
tinction between perfect and imperfect frustration in the 
situation (B) of inequivalent layers. 



D. Zero-field phase transitions 

For perfect frustration, upon increasing the ratios J'/J 
(e.g. by applying pressure) either (pA or may condense 
first. If the difference between the masses, to^b — m^A, 



Similar to the zero-field transition, the field-driven 
transitions at Hd and Hc2 are split in the presence of 
perfect frustration due to the absence of a proximity ef- 
fect. The primary A ordering transition is in the dilute 
Bose gas (i.e. BEC) universality class with z = 2. At 
Hci, the dimensional crossover between 3d and 2d criti- 




FIG. 15: a) Schematic overall temperature-field phase dia- 
grams for the coupled-dimer model with inequivalent layers 
and perfect frustration, for different relations between Ja.b, 
J'a b- Solid (dashed) lines show the phase transitions on the 
A (B) subsystem, b) Schematic zero-temperature magnetiza- 
tion curves for the three scenarios in a). For perfect frustra- 
tion, the slope will change discontinuously at the transitions. 
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FIG. 16: a) Triplon energy gaps and b) uniform (solid) and 
staggered (dashed) magnetization as in Fig. 9, but for inequiv- 
alent layers with perfect frustration - the splitting of the phase 
transitions is clearly visible. Parameter values are: Ja ~ 4.27 
meV, Jb = 5.04 meV, J'a = J'b^ 0.5 meV, -J2z = Jz = 0.1 
meV. The field is displayed in units of the average J = 4.66 
meV. The thin lines in b) show the individual contributions 
of the A and B subsystems (i.e. even and odd layers) to the 
magnetizations. The presence of two Goldstone modes in the 
intermediate field range is an artifact of the harmonic approx- 
imation: inclusion of interaction effects would split the two 
modes. 



cal behavior is again set by the scale 77^ . (As discussed in 
Sec. VD, this scale vanishes at Hc2, if the Hamiltonian 
has no bare J^z ) 

Once the A subsystem is ordered, the ui term intro- 
duces an easy-axis anisotropy for the remaining <f>B U(l) 
degrees of freedom (irrespective of the sign of ui ) . Hence, 
the secondary transition at H'^-^ is in the Ising universality 
class with z = 1. (The change from z = 2 to z = 1 can be 
understood, e.g., in terms of triplon modes: the easy-axis 
term couples and r_, and the resulting 2x2 matrix 
diagonalization results in a linearly dispersing mode at 
criticality. The easy-axis term is neglected in the har- 
monic bond-operator calculation.) The crossover with 
increasing energy is now from 3d Ising to 2d Ising to 2d 
BEG, with crossover scales as above. 

This discussion implies that, for perfect frustration, 
there will be a field range Hd < H < H'^^ where the 
B layers are still disordered. In this regime, the B layer 
magnetization is zero (exponentially small) at T = (low 
T), respectively. This will only change at the secondary 
transition which occurs where the B gap as function of 
the field reaches zero, leading to a kink e.g. in the total 
magnetization, see Figs. 15 and 16. Note that this kink 



FIG. 17: a) Triplon energy gaps and b) uniform (solid) and 
staggered (dashed) magnetization as in Fig. 16, but now for 
inequivalent layers with imperfect frustration - here the sec- 
ondary phase transitions are smeared out due to the proximity 
effect. Parameter values are: J a = 4.27 meV, Jb = 5.04 meV, 
J'a^ J'b= 0.5 meV, J^i = 0.05 meV, = 0.15 meV - the 
latter two values are the couplings along the two inequivalent 
Jz diagonals, see Fig. lb. 

will also be present at finite T, due to the power-law 
onset of the secondary order parameter at the transition. 
The distance between the two transitions is essentially 
given by the difference in the zero-field gaps of the two 
modes 0a, (t>B\ small corrections arise from the density 
interactions between the two condensates. 

As above, for imperfect frustration, the proximity ef- 
fect will smear out the secondary transition, and with it 
the magnetization kink. The energy scale for the dimen- 
sional crossover is determined by the bare dispersion of 
the lowest triplon and scales as J^. 

Magnetization curves obtained using the bond- 
operator method are shown in Figs. 16 and 17, for per- 
fect and imperfect frustration, respectively. For imper- 
fect frustration, the magnetization tob of the B subsys- 
tem is now finite above Hd , but depends in a non-linear 
fashion on the A magnetization mA (or the field). 



F. Larger unit cells 

For unit cells containing more than two dimers, the 
number of triplon modes increases accordingly, but many 
qualitative features of our analysis remain valid, (i) 
There will be at maximum two distinct condensates, not 
coupled by a proximity effect in the case of perfect frus- 
tration. Hence, the number of separate sharp ordering 
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transitions will not exceed two, although the number of 
mode gaps is larger, the reason being a proximity effect 
due to unfrustrated couplings, (ii) The bandwidth of the 
vertical dispersion for N inequivalent layers will in prin- 
ciple be determined by iV-th order perturbation theory; 
depending on the relation between the Jz and the mode 
energy differences Awq, bandwidths between cx Jz and 
oc are possible. 



VII. APPLICATION TO BACUSI2O6 

After having collected information about the different 
scenarios, we finally come to a more detailed discussion 
of the behavior of BaCuSi206. First we try to summarize 
important experimental results: (i) The phase boundary 
near Hd = 23.2 T follows the 2d power law with shift ex- 
ponent -ip = 1 down to 30 mK.^®''*'* (ii) Zero-field neutron 
scattering shows multiple modes, implying inequivalent 
dimers.^^ Together with the structural distortion below 
100 K this suggests inequivalent layers. Further, the ver- 
tical dispersion of the modes seems tiny (< 0.1 meV) 
even at wavevectors away from the frustrated point, (iii) 
NMR experiments show the presence of two inequivalent 
Cu sites above Hd, with significantly different, but non- 
zero magnetizations.^^ (iv) The magnetization curve is 
almost linear above Hd-^^ (v) Small spin anisotropics 
arise from magnetic dipolar interactions, but those are of 
order 10 mK only."*^ 

Points (ii) and (iii) rule out a perfect bet lattice struc- 
ture with a single-dimer unit cell. Sec. V. In the following 
we therefore assume the inequivalence of even and odd 
layers, i.e., our scenario (B) of Sec. VI. Further, we will 
ignore the possible tendency towards incommensuration, 
i.e., we stick to layers with (tt, tt) ordering. We are aware 
that this may not fully account for all details relevant 
to the material. (Note that a full characterization of the 
low-temperature lattice structure is still open.) 

Within the scenario of inequivalent layers, perfect frus- 
tration implies two distinct condensates without proxim- 
ity effect. This would invariably lead to two transitions as 
function of field (even at finite T) , with a kink in the mag- 
netization curve. The second transition field can be esti- 
mated from the zero-field gaps^^ to be around H'^^ « 27 
T. In the field range Hd < H < H'd the magnetization 
in the second (non-ordered) subsystem should be expo- 
nentially small. In contrast, the NMR data^'^ indicate 
a sizeable magnetization here and hence appear incon- 
sistent with the perfect-frustration scenario. Further, to 
our knowledge no secondary thermodynamic transition 
has been observed near 27 T. Together, this appears to 
rule out perfect frustration. 

We are left to consider inequivalent layers with im- 
perfect frustration. Imperfect frustration induces a fi- 
nite vertical bandwidth for all in-plane wavevectors, and 
it causes proximity effect, thereby smearing out the sec- 
ondary ordering transition. The finite vertical bandwidth 
at 5|| = (tt, tt) is the relevant dimensional crossover scale 
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FIG. 18: Zero-field triplon dispersions for inequivalent lay- 
ers with partially frustrated and modulated vertical coupling, 
Fig. Ic. Parameter values are: Ja = 4.27 meV, Jb = 5.04 
meV, J'j^ = J'g = 0.5 meV, J^ai = 0.09 meV, J^a2 = 0.11 
meV, Jzsi = 0.01 meV, JzS2 = 0.19 meV. The bandwidth 
along (tt, tt, q^) is reduced by a factor of 3 as compared to the 
situation without modulation (Fig. 11) - and can be made ar- 
bitrarily small by reducing \ JzAi — JzA2\, but the magnetiza- 
tion data in the present case are essentially undistinguishable 
from the ones in Fig. 17 above. 

below which 3d behavior is observed at the phase transi- 
tion. From (i) , this scale is known to be less than 30 mK - 
this significantly constrains the unfrustrated component, 
Jzi — Jz2, of the inter-layer coupling (Sec. VIA): it has 
to be smaller than 0.035 meV (assuming this coupling to 
be vertically unmodulated; note that the square of this 
value is proportional to the vertical bandwidth). Tak- 
ing this value, we have calculated magnetization curves 
which only display tiny deviations from the case of per- 
fect frustration. Fig. 16, i.e., a pronounced magnetization 
kink at a second "transition" field ff^j^, and essentially 
zero magnetization in the B subsystem below this field. 
Apparently, the proximity effect is too small to be in 
agreement with both the NMR and the magnetization 
data. 

Thus, a 3d scale smaller than 30 mK is not easily com- 
patible with the large proximity effect. To reconcile these 
facts, we propose the inter-layer coupling to be both frus- 
trated and vertically modulated, see Fig. Ic. In such a 
situation, the unfrustrated components of the inter-layer 
coupling can alternate between being weak and strong, 
rendering a large proximity effect compatible with a tiny 
3d dispersion. To prove the case, we have calculated 
dispersions and magnetization data in a situation with 
JzAi = 0.11 meV, JzA2 = 0.09 meV, JzBi = 0.19 meV, 
JzB2 = 0.01 meV, i.e., the unfrustrated component of the 
inter-layer alternates strongly, while the frustrated com- 
ponent is unmodulated. With this choice, we obtain mag- 
netization data with a sizeable proximity effect very sim- 
ilar to Fig. 17, but the vertical dispersion along (tt, tt, qz) 
[which is given by {JzAi - JzA2){JzBi - JzB2)/^t^q] is 
strongly reduced, Fig. 18. By reducing \ JzAi — JzA2\ fur- 
ther, this 3d scale gets arbitrarily small, without signifi- 
cant changes in the magnetizations. (Note that the pre- 
cise shape of the magnetization curve will be infiuenced 
by logarithmic corrections not captured by our harmonic 
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approximation, but the presence or absence of a rather 
sharp kink will be unaffected. As stated above, we never 
find a field-independent ratio of the two magnetizations, 
m A I TUB, and attribute the experimental observation^'^ of 
rriA/mB « 5 to a coincidence related to proximity effect 
and logarithmic corrections.) 

In summary, within the scenario of modulated verti- 
cal couplings we can easily get consistency between two 
key observations: a tiny 3d scale and a large proximity 
effect. We note that a modulation of the unfrustrated 
part of the coupling is required for both phenomena; the 
largest unfrustrated coupling is of order 0.2 mcV (Fig. 18) 
and hence relevant for the magnetism not only at low- 
est T. There are more free parameters in the proposed 
model: an additional modulation of the frustrated part 
will mainly reduce the bandwidth along (0, 0, qz) - here 
future neutron scattering experiments will help to put 
bounds on these couplings. In the absence of further 
experimental information on the 3d dispersion scales we 
refrain from a more detailed fitting of the experimental 
data. 

Let us finally comment on the transition temperature 
Tn. In a general quasi-2d XY magnet with some (unfrus- 
trated) inter-layer coupling J^, Tn will sensitively depend 
on Jz'. For Jz — > 0, Tn is given by the Kosterlitz-Thouless 
temperature Tkt, but a finite Jz will induce a strong (log- 
arithmic) enhancement to that, due to the exponentially 
large correlation length above the 2d Kosterlitz-Thouless 
transition. For a partially frustrated inter-layer coupling, 
both the unfrustrated and frustrated parts will contribute 
to Tn, but the detailed behavior may be complicated (see 
the discussion in Sec. VE), and a quantitative calcula- 
tion of Tn is beyond the scope of the paper. Applied to 
BaCuSi206, we believe that the rather small value of the 
unfrustrated 3d dispersion scale (< 30 mK) may well be 
compatible with the observed Tn (which is roughly twice 
as large as the expected Tkt), as the frustrated part of 
the inter-layer coupling is likely significantly larger (< 0.5 
K) and enhances Tn as well. 



VIII. CONCLUSIONS 

We have studied quantum dimer magnets on layered 
lattices of the bet type. Guided by experimental results 
on BaCuSi206, we have discussed two distinct routes 
towards magnetism with reduced dimensionality: (A) 
an ideal bet lattice with perfect frustration, and (B) 
a distorted lattice with inequivalent layers and possi- 
bly broken frustration. While the situations with per- 
fect frustration arc theoretically fascinating, the mate- 
rial BaCuSi206 likely falls into class (B) without per- 
fect frustration. There, dimensional reduction primarily 
arises from a layer "mismatch" : vertical magnon hopping 
proceeds via a second-order process of strength J^/Awg 
(where Atjq is the magnon energy difference between 
even and odd layers). In addition, we propose that a 
further reduction of the 3d energy scale originates from a 



modulation of Jz along the c axis, with the unfrustrated 
part being small on every second link. This scenario al- 
lows to reconcile the presence of inequivalent layers^"^'^^ 
and the tiny 3d energy scale^* with the rather large mag- 
netic proximity effect. 

For the case of perfect frustration, a number of theo- 
retical challenges are left for future work: Those include 
(i) a more detailed study of the dimensional crossover be- 
havior in the quantum critical regime, and (ii) the fate of 
the order-from-disorder mechanism at elevated tempera- 
tures and the nature of the finite-temperature transition, 
including a reliable estimate of the Neel temperature. 
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APPENDIX A: BOND OPERATORS BEYOND 
THE HARMONIC APPROXIMATION 

In this appendix we present details of the paramagnetic 
bond-operator approach introduced in Sec. IV B - this 
is a generalization of the formalism of Ref. 28 to finite 
fields. An additional generalization to finite temperature 
is possible as well, but we shall not follow that here. 

The Hamiltonian is given by W = W2 + 'H4 + Hu, with 
H2 and Hu given in (15) and (20), respectively. The 
quartic term 7^4 = Ti^i + H4Z is 



Ha 



J' 

— hi{in;jr 



T^iz = h4{in; i -|- A, n -|- 1), 



(Al) 



iAn 



with dimer site indices {in) as in Eq. (1), J4Z — J}, 
Jf^ -I- J_J^ -I- J|^, and /14 given by 



hi{i;j) 



ft t I tt t 



(A2) 



-ft (ff^+ f 



h.c. 



In principle, there is also a cubic term. 



iAn 



j+ j- 

-^h3{in;i + A,n+l) + ^h3{i + A,n+l;in) 

(A3) 
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with jf^ = Jl^-Jf±{Jl''-Jf) and 



a = +): 



iAn 

[(t,+ +4-)(4oV-*]+*.-o) + H)- 



+ 



We neglect it as its inclusion would not qualitatively af- 
fect our conclusions: Because of the geometric frustra- 
tion, self-energy contributions deriving from Ti^z which 
are cx Jf^ (Fig- 4d) give no gz-dependent correction to 
the triplon dispersion for gj| = (7r,7r). 

The effect of the local and non-retarded interaction 
Hu can be captured by a ladder summation, leading to 
a four-point vertex 



-y 



P fe— p 



(A5) 



where N is the number of lattice sites. Here, all anoma- 
lous scattering vertices have been neglected, which is jus- 
tified if the number of triplet bosons is small (Brueckner 
approximation). (Note that the Bogoliubov coefficients 
u^v do not depend on the Zeeman index, see below.) 

To leading order, the normal self energy from Jiu is 
given by the sum of Hartree and Fock diagrams (Fig. 4a): 

I]^(fc,w) = I]„a(^,w) + X]^«/3(^''^)' 



T,a)3{k,0j) = —^V^i3Tap,ap{k + q,UJ-UJ^i3).{A6) 



For the bet lattice, this self-energy will induce a triplon 

dispersion along (7r,7r, (72) of order J|, with the leading 
process being the one in Fig. 3b2 (the scattering ampli- 
tude from Hu is of order unity). 

In addition to the hard-core interaction Hu, we take 
into account the quartic t terms in H4 in a mean-field 
approximation, Fig. 4b. This gives additional self ener- 
gies: 



with 7g^| and 7^^ defined in Eq. (16). While these mean- 
field terms do not contribute to the dispersion along 

(tt, TT, (jx), higher-order processes in Haz do so. Remark- 
ably, the second-order self energy. Fig. 4c, turns out to 
give the leading contribution to the vertical dispersion, 
scaling with jf (J^ from the vertices and for vertical 
hopping on one internal line). Therefore we include a 
contribution to S"(fc), which is of second order in H^z, 
in addition to the Hartree-like terms above. The evalu- 
ation of the diagram in Fig. 4c gives (here displayed for 



E---(fc,a;) = Jl^,Y.^UAk,,^+,m 



pq 



+2%-,-.,(o, 0, -) + +, +) + 4..^(+, -) 

(A8) 



with 

fkpqj^'f^'^y) 



2 9 9 

k+p-q <1 P 



(A9) 



The normal self energies S'^ and S" add up, S"(fc) = 
S"(A?) + YFik^LO = 0). With the self energies we can 
define renormalized coefficients A^:, with 

Aj: = Aj: + S"(fc), Bj: = Bj: + S'^(fc). (AlO) 

Solving the coupled Dyson equations for the normal 
and anomalous Green's func;tions yields the renormalized 
spectrum for the quasiparticles, 



Or: = ZrjAi- Bi-aH, 

ka ky k k ' 

with a quasiparticle weight 

aE^(fc,w) 



dto 



(All) 



(A12) 



w=0 



and the renormalized Bogoliubov coefficients: 



±- 



1 



.1, 



UkVk = - 



Bk 



(A13) 



A fully self- consistent solution of the problem requires 
the substitutions 



^ka 



'ZkU^k' 



ZkVk 



(A14) 



in the equations (A5,A6,A7,A9). 

The reader may wonder about a possible a (spin) de- 
pendence of the self energies and the Bogoliubov coeffi- 
cients. Indeed, the naive perturbative approach sketched 
here results in such a spin dependence, which conse- 
quently also yields a rcnormalization of the magnetic field 
(or the g factor). Physically, such a result is incorrect, 
as the z component of the total spin is conserved. This 
conservation law can be formally implemented by a corre- 
sponding Ward identity, resulting in appropriate vertex 
corrections. Here, a consistent result can be obtained 
by implementing a finite field through the replacement 
CO ^ Lo + aH everywhere. This implies that the self- 
energies need to be evaluated at energy aH instead of 
zero. Doing this, the self energies and Bogoliubov co- 
efficients remain a independent."'^ (Practically, only the 
a = self energy needs to be evaluated.) 
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APPENDIX B: SCATTERING AMPLITUDE FOR 
BETHE-SALPETER EQUATION 



Altogether, the scattering amphtude in the singlet {S = 
0) channel is 



We consider a biquadratic coupling between adjacent 
layers, i.e. a simplified version of Eq. (4): 



(Bl) 



rewrite it in terms of bond operators and extract from 
the quartic terms the following contribution to the bare 
scattering amplitude 



» j^coll 



J. 



coll 



(B2) 



where we employ the notation of Ref. 37. This adds to 
the result deriving from Hi and Hu, 

+2(^a5^/37 - 5al3^-i5){J' lp-q,\\ + JAzlp-q,z) 

+U {6as5p^ + SaS^p^). (B3) 



-4J'(7p+g,|| + 7p-9-,||) - 4:J4z{7p+q,z + lp-q,z) 

+2U + Jcoii(27p-+g-,^ + '2^p-q,z + 7qJ- (B4) 



Interestingly, the quartic scattering terms arising from 
Jiz and Jcoii are of essentially identical form (although 
the Hamiltonians are different). Consequently, both Haz 
and WcoU cause similar triplon bound states. 
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